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The atomistic structure and energetics of the E13 (10l4) 
symmetrical tilt grain boundary in a-Al2 03 are studied by 
first-principles calculations based on the 
local-density-functional theory with a mixed-basis 
pseudopotential method. Three configurations, stable with 
respect to intergranular cleavage, are identified: one 
Al-terminated glide-mirror twin boundary, and two 
0-terminated twin boundaries, with glide-mirror and 
two-fold screw-rotation symmetries, respectively. Their 
relative energetics as a function of axial grain separation are 
described, and the local electronic structure and bonding are 
analysed. The Al-terminated variant is predicted to be the 
most stable one, confirming previous empirical calculations, 
but in contrast with high-resolution transmission electron 
microscopy observations on high-purity diffusion-bonded 
bicrystals, which resulted in an O-terminated structure. An 
explanation of this discrepancy is proposed, based on the 
different relative energetics of the internal interfaces with 
respect to the free surfaces. 
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I. INTRODUCTION 

Alumina (q;-A1203) is widely used in polycrystalline 
forms. Small amount of various dopant elements (typi- 
cally Mg, Ca, Si, and Ti) are added in the synthesis of 
alumina polycrystals, either to promote sintering, or to 
obtain empirically optimized microstructures (grain size, 
and shape distributions) for technological applications.EI 
Due to the very small solubility limit in bulk for npany 
impurities (typical values are 10 ppnufor yttrium,0 30 
ppm for calcium,El 300 ppiaa for siliconQ), they segregate 
to grain boundaries (GB).a The segregation of different 
species at grain boundaries drastically modifies the prop- 
erties of alumina, in either beneficial or undesired man- 
ner. For instance, the isovalent cations Y and La reduce 
the mechanical creep rate poi polycrystalline alumina by 
2—3 orders of magnitiideJj'El Si and Ca promote exag- 
gerated grain growth,BLl and Ca may affect the fracture 
properties decreasing the strength and toughness of the 
materiaLml From a fundamental aspect of science, a key 
issue in alumina systems is therefore the microscopic un- 
derstanding of internal interfaces and of grain-boundary 
impurity segregation. 

This study is part of an extended project that 
aims to correlate the macroscopic structural proper- 



ties of alumina to its microscopic structure. The 
project combines empirical atomistic modelling and 
first-principles electronic-structure theory with experi- 
mental high-resolution transmission electron microscopy 
(HRTEM) and with analytical electron microscopy 
(AEM), in particular spatially resolved electron energy- 
loss near-edge spectroscopy (ELNES). In this framework 
a set of four grain boundaries were selected for a fun- 
damental theoretical study of atomistic and electronic 
structures: the basal S3 (0001), the rhombohedral S7 
(1012), the E13 (1014), and the prismatic E3 (1010) twin 
boundaries. This set provides a couple of mutually or- 
thogonal twin interfaces, basal and prismatic, and a cou- 
ple of nearly orthogonal ones, the E7 and the S13. 

The basal and rhombohedral twin interfaces have re- 
cently been studied in detail by means of ab-initio local- 
density-functional tiieacy (LDFT) and empirical shell- 
model calculations.EJ"Ej The theoretical analysis of the 
interface structure for the rhombohedral twin has been 
confirmed quaptitatively by both HRTEM and ELNES 
experimentsjl3 validating the accuracy and predictive 
power of the theoretical approach, which is employed as 
well in the present work dealing with the S13 twin in- 
terface. Results forpJjhe prismatic twin interface will be 
reported separatelyEil. 

Several studies indicated that the segregation of im- 
purities like Ca or Y to GB is more favourable for gen- 
eral, disordered interfaces than for special, opes with high 
structural order and energetic stability.E3il3 This is com- 
monly rationalized by the large mismatch in ionic radii 
between the segregated species (r,=0.94 A for Ca2+ and 
ri=0.88 A for Y^+, Ref. |l^) and the host-lattice ele- 
ment Al^+ (ri=0.45 A, Ref. |l^). General boundaries 
provide more excess volume at the interface, and the 
oversized Ca and Y cations fit more easily here than at 
special GB, where the interfacial atomic environment is 
highly ordered and compact. As an example, the impu- 
rity level in the nearly bulk-like rhombohedral and pris- 
matic twin baundaries is below the detection limit of 0.3 
atoms/nni'^.ll3 There are indications that the basal plane 
may have an atomic envicoppient which is slightly more 
favourable to segregation.E3~c2l Compelling experimental 
evidence for segregation to a twin boundary waLS,observed 
for the case of the E13 (10l4) twin boundaryH 

Among the four selected GB, the interfacial structure 
of the S13 boundary is the most different from the bulk 
atomic environment of q;-A12 03. Hence it is anticipated 
to have a rather high interface energy, and it seems to 
be best suited and most promising to address the issue 
of segregation of cation impurities to GB. Its interfa- 
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cial structure appears sufficiently disordered and open 
to accept likely segregated impurities, thus being repre- 
sentative for the atomic environment in a general GB. 
At the same time it is sufficiently ordered, with short 
lateral periodicity along the interface, that it can be de- 
scribed with some tens of atoms. Hence it is accessible via 
ffist-principles LDFT calculations. In the present work 
we consider the S13 twin boundary as a representative 
case to study the effect of segregation on the mechanical 
and electronic properties of alumina interfaces. We de- 
scribe the properties of the pure S13 twin interface in this 
paper, and the modifications induced by the segregated 
atoms in a forthcoming one. 

Theoretical investigations preceding ours of the same 
E13 twin have been undertaken by Finnis and coworkers. 
To our best knowledge, results have beea Jiiiblished only 
for empirical shell- model calculations ;l2I EJ but not for 
ab-initio LDFT calculations. Ej 

This paper is arranged as follows: A short review of 
previous studies on the S13 twin boundary is presented 
in Section |l[ A brief desciption of the employed compu- 
tational methods is given in Sec. 111. The results of our 



investigation, namely the mechanical stability, the rela- 
tive energetics, and the electronic structure of different 
metastable structural variants for the E13 interface are 
reported and discussed in Sec.JfV[ Final remarks and a 
summary are made in Section |V[ 



II. PREVIOUS STUDIES ON THE E13 TWIN 
BOUNDARY 

Figure 1^ shows the crystal structure of a-Al2 03 (corun- 
dum, space group i?3c), projected on the (1210) plane: 
light and dark circles represent oxygen and aluminium 
ions, respectively. The structure can be described as a 
nearly ideally close packed hexagonal oxygen sublattice 
with 2/3 of the available octahedral interstices filled by 
Al cations. The remaining octahedral sites are empty, 
and the partial filling results in small internal atomic 
displacements in both sublattices, which are visible in 
Figure 0. The E13 (1014) twin interface is formed by 
two grains with a relative misorientation given by the 
planes (10i4)||(10l4) and the directions [1210]||[12l0]. 

High-resolution transmission electron micrographs of 
high-puriti' diffusion-bonded bicrystals containing this 
interfaca23 revealed two distinct boundary structures. 
The images differed by the sequence of intensity maxima 
across the boundary: one showed a continuous (C) and 
the other a broken (B) pattern of spots across the bound- 
ary. The C interface was found in the bicrystal bonded 
without further treatment, besides sufficient cleaning of 
the free surfaces, while the B and C interfaces coexisted 
when a metallic Al layer was deposited on the free sur- 
faces after the cleaning and before the diffusion-bonding 
process. 

Atoraisfe 
tentialsE3S 



modelling based on empirical classical po- 
also identified two boundary structures, in 



which the two grains were either Al or O terminated. 
They predicted the former to be more stable than the 
latter. The calculated interface energies for the two re- 
laxed boundaries were 1.7 J m~^ (Al-terminated) and 2.1 
J m^^ (0-terminated), respectively. 

The details of these interface structures were not di- 
rectly accessible from the HRTEM images alone. The 
combination of atomistic modelling and image simula- 
tions allowed an assignment of the B and C patterns ta 
cation and anion terminated interfaces, respectively.tj 
The corresponding terminating planes in the perfect 
corundum structure of Figure |l| are marked by the dashed 
lines labelled B and C. 

Hoche and Riihle studied the segregation of calcium 
to the E13 boundary and the resulting structural modi- 
fications induced on the B and C structures. u A layer of 
metallic calcium, 3 nm thick, was deposited on one of the 
(1014) free surfaces after cleaning and before diffusion- 
bonding the bicrystal. After the bonding proceduKe, 
2.4±0.9 Ca ions per nm^ were detected at the interface,E£l 
therefore proving that Ca segregates to this grain bound- 
ary. Again, as in the Al-doped case, two distinct struc- 
tures coexisted: one variant showed little differences to 
the C structure. The other was qualitatively similar to 
the B structure but with significant quantitative differ- 
ences. In this case the comparison of the HRTEM with 
theoretical calculations was not possible because atom- 
istic calculations of cation segregation at alumina grain 
boundaries were not available. 

Such results raise fundamental questions on the for- 
mation mechanism and kinetics of interfaces. In the 
specific case of the E13 twin boundary, the most com- 
mon structure observed in diffusion-bonded bicrystals 
(0-terminated) is apparently not the one with the min- 
imum energy predicted by the theoretical calculations 
(Al-terminated). Hence, is the E13 twin boundary a case 
in which empirical atomistic model potentials reach their 
limit of predictive power? This concern is assessed by em- 
ploying the non-empirical LDFT approach for the theo- 
retical investigation. Or is the E13 twin boundary rather 
a case where diffusion-bonding of the lowest-energy sur- 
faces does not immediately result in the lowest-energy in- 
terface? This possibility is investigated by mapping theo- 
retically the relative energetics along configuration paths 
for both bonding the most stable surfaces together and 
pulling the most stable interface apart. The paper pro- 
vides a precise description of the interfacial atomic envi- 
ronment and electronic structure of the E13 twin bound- 
ary, and it proposes an explanation of the difference be- 
tween experimental HRTEM observations and theoretical 
predictions. 



III. COMPUTATIONAL METHODS 

In this work we used the same computational apjjtoach 
as previously for the rhombohedral twin boundary,EiJ the 
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basal twin boundaryjlj and basal-plane stacking faults,EZI 
namely atomistic modelling with empirical interatomic 
potentials and ab-initio LDFT calculations. 

The empirical potentials which were used in the atom- 
istic modelling of a-Al2 03, were taken from Lewis and 
CatlowEa. They incorporate the long-range Coulomb in- 
teraction between all Al cations and O anions, the short- 
range ion-overlap repulsion, the weak Van-der-Waals at- 
traction between the anions, and a shell-model dipole po- 
larizability of the anions. The cations were treated as non 
polarizable and mutually interacting via Coulomb repul- 
sion only. The advantage of the empirical potentials with 
respect to LDFT is their superior computational simplic- 
ity, the disadvantage is their limited reliability for pre- 
dictions of material-specific properties. The comparison 
of the two approachestilEZi for basal-plane twin-boundary 
structures and stacking faults in a-Al203 showed that the 
empirical potentials can provide quantitatively not very 
accurate interface energies, but qualitatively correct en- 
ergy hyperfaces in structural parameter spaces. Atom- 
istic arrangements in metastable structures obtained by 
both approaches were found to be very similar, making 
the empirical potentials very useful for identifying rele- 
vant initial grain-boundary configurations, that can then 
be structurally optimized with the LDFT. 

The LDFT calculations were made by means of the 
mixed-basis pseudopotential technique for the computa- 
tion of electronic structures, total energies, and forces 
in crystalline materials. The interfaces were described 
by supercells with periodic boundary conditions. The 
computational details, namely the norm-conserving ionic 
pseudopotentials for Al and O, the mixed basis for the 
representation of the valence-band and conduction-band 
Bloch states in q;-A1203, the discrete Brillouin-zone in- 
tegration meshes etc. were the same as, or eouiiialent 
to the ones used in our previous investigations.lly~t2l For 
the optimization of the interface structures, excess ex- 
pansions of the interfaces were obtained via total-energy 
minimization with respect to the axial supercell length. 
Simultaneously, to account for lateral grain translations 
as well, all atomic positions in the supercells were stati- 
cally relaxed according to the acting forces. 

Numerous published studies of various materials' 
structures and properties from the Stuttgart groups and 
others involved in deV|ek)|ptnents of the mixed-basis pseu- 
dopotential approachnjO have documented the accu- 
racy, efficiency and predictive power of this density- 
functional technique. Specifically for q;-A1203, bulk 
structural parameters of corundum were compared to re- 
sults from other LDFT methods and from experiment 
in Table 1 of Ref. and the predictive power of the 
method for grain-boundary structures was demonstrated 
in Ref. 1^. 

The unit cells used to model the E13 twin interface 
were constructed by cutting the corundum structure of 
a-Al203 with a (1014) plane and rotating one half of 
the crystal with respect to the other by 180°, around the 
interface-plane normal vector nGB = [5054] (symmetrical 



twist). Note that in hexagonal coordinates the [5054] 
direction is orthogonal to the (1014) plane. The rota- 
tion axis ncB contains a center of inversion symmetry 
for the original corundum structure, therefore the result- 
ing bicrystal is mirror-symmetric. 

This procedure was applied to the cutting planes B and 
C marked by the dashed lines in Figure |l|, leading to the 
starting configurations for two stoichiometric variants of 
the grain boundary. The computational cell was then 
obtained by cutting a slab centered on the grain bound- 
ary, and repeating it with periodic boundary conditions. 
The slab-edges identifying the shape of the orthorhombic 
superceh are ei || [2021], 63 || [1210], and 63 || [5054]. 

The size of the slab along 63 was selected so that it 
contained two equivalent interfaces, one in the middle, 
already described, and one at the edges, arising as a con- 
sequence of the periodic boundary conditions. A care- 
ful choice of the supercell size along 63 was necessary in 
order to minimize the spurious interaction between the 
two equivalent boundaries. The convergence of the re- 
sults against the interface-interface distance was checked 
by considering 3 supercells, with 10, 15, and 20 atomic 
layers between the boundaries, corresponding to 40-, 60- 
, and 80-atom supercells. The results for the small cell 
were already converged in the case of the cation termi- 
nated interface while the 60-atom cell was necessary for 
the anion terminated one. We will further discuss this 
issue in the following. 

IV. RESULTS AND DISCUSSION 
A. Lateral translation state 

The twin-boundary equilibrium configurations were 
found with increasing degrees of accuracy, by combin- 
ing empirical atomistic and first-principles calculations. 
Previous studies showed that the empirical shell-model 
potential provides the correct qualitative description of 
the energetics and structure of twin boundaries. t3'tll This 
makes the shell-model suitable for a first exploration of 
possible metastable configurations. j— . 

Using this model and the computer program GULP,Lj 
we calculated the bicrystal energetics as a function of 
the relative displacement of one grain with respect to 
the other along the interface plane. The relative dis- 
placement is defined by the lateral translation vector 
TiGi + 1262, where Ti and T2 are fractional coordinates 
with respect to the lattice vectors ei and 62. The inter- 
facial symmetry confines the non-equivalent part of the 
energy surface to a quarter of the lateral repeat unit cell, 
limited by the interval between and 0.5 for both Ti and 
T2. 

The starting boundary structures, corresponding to 
(Ti, r2)=(0,0), were the mirror-symmetric interfaces ob- 
tained as described in the previous section. It is clear 
that these structures are highly unstable because equal 



3 



charges on opposite sides of the boundary face each other 
at very close distances. This is the reason why atomic 
relaxation was not considered in this preliminary analy- 
sis of the translation state. The shell-model potentials 
are not suitable to describe such highly non-equilibrium 
structures and would predict large unrealistic relaxations. 
Only the metastable structures identified with the static 
shell-model calculations were relaxed subsequently, ac- 
cording to both empirical shell-model and first-principles 
forces on the individual atoms. 



1. Cation terminated interface 

The cation terminated interface was obtained by cut- 
ting the corundum structure along line B in Figure land 
building the supercell as described above. The energy 
surface, governing the inter-granular lateral translation 
state, as a function of (ri,T2) is shown in Figure ||. 

The starting configuration (Ti, T2)=(0,0) is expected 
to be unstable because it faces equal charges very close 
to each other, on opposite sides of the boundary. The 
shell-model predicts the highest energy for this config- 
uration. The total energy is minimized by the lateral 
translation state (Ti, r2)=(0.25,0.5), resulting in an Al 
terminated structure with glide-mirror symmetry with 
respect to (1014). This configuration, which we will in- 
dicate as G(A1), maximizes the inter-granular distance 
between equal charges and minimizes that one between 
opposite charges. Alternating triplets of oxygen atoms 
at the interface (0(l)-0(2)-0(3) and 0(l)*-0(2)*-0(3)* 
in Figure ^) characterize this boundary. 

The structure was then relaxed according to the forces 
from the shell model as well as from the LDFT. The 
results were qualitatively equivalent: the stable transla- 
tion state was the same, and internal relaxations involved 
only ions in the neighbourhood of the interface, namely 
the atomic layers from 0(2) to 0(2).* The LDFT relaxed 
structure is shown in Figure ^. It has the same features 
as the cation-terminated interface identified in Ref. |2^: 
an HRTEM setting for this interface showing the holes 
as intensity maxima would produce a broken pattern of 
white spots across the boundary. 

In order to check the adequacy of the supercell size 
perpendicular to the interface, the relaxation procedure 
was made for 40- and 60-atom supercells. The small dif- 
ference between the interface energies (< 4%) of the two 
supercells showed that convergence was already achieved 
with the small cell. 



2. Anion terminated interfaces 

A similar procedure was carried out for the anion ter- 
minated interface, which was obtained by cutting the sin- 
gle crystal along line C in Figure 0. Note that in this case 
the interface plane coincides with an atomic oxygen layer. 



The total energy as a function of the lateral translation 
state {Ti,T2) is shown in Figure | and predicts three 
metastable configurations. 

The structure with the lowest energy has a transla- 
tion state of (Ti, T2)=(0.25,0), and corresponds to the 
Al-terminated boundary already described above. This 
is an important feature of this boundary: the same Al- 
terminated interface may be obtained either by cutting 
the single crystal along B (Figure |l|), rotating one grain 
around ncB by 180° and imposing the lateral translation 
state (Ti, T2)=(0.25,0.5), or by cutting along C, rotating 
one grain around nQB by 180° and imposing the lateral 
translation state (Ti, T2)=(0.25,0). In the latter case the 
initial boundary was on an O layer, but the glide-mirror 
symmetry operation reconstructs a bulk-like environment 
near the atomic layer at C, and shifts the boundary by 
two atomic layers from O termination (C) to Al termi- 
nation (B). 

The structure with the second lowest energy is defined 
by the lateral translation state (Ti, r2)=(0.5,0.5). It has 
an oxygen layer on the boundary plane, and the grains 
are related by a glide-mirror symmetry with respect to 
(10l4). We will denote it by G(0). Similarly to the Al- 
terminated case, both the LDFT and shell-model atomic 
relaxations predict the metastability of this translation 
state: the atomic relaxation involves mainly the first 
three interfacial layers. The relaxed supercell is shown 
in Figure ||. Note the characteristic 'V shape formed by 
the interfacial O atoms (O(l)-O(0)-O(l)* in Figure |). A 
HRTEM micrograph of this structure showing the holes 
as white spots would produce a continuous pattern across 
the boundary. 

Finally, the metastable structure with the highest en- 
ergy is obtained by the translation state (Ti, r2)=(0,0.5), 
leading to an 0-terminated configuration with a twofold 
screw-rotation symmetry around [2021]. It is denoted as 
S(0). The relaxed supercell is shown in Figure ^. Both 
the LDFT and shell-model atomic relaxations show that 
this translation state is not stable. The atomic forces 
drive the bicrystal towards (Ti, T2)=(0.12,0.5). This is a 
new metastable interface structure that was not noticed 
in the previous studies of this GB. Holes would produce 
a broken pattern of spots, intermediate between the G(0) 
and G(A1) ones. 

Also for G(0), the adequacy of the axial supercell 
size was checked by comparing the relaxation patterns 
obtained with 40-, 60-, and 80-atom supercells. The 
LDFT relaxation of the 40-atom cell involved all the 
atoms hence leaving virtually no bulk-like volume re- 
mote from the boundary, while the shell-model relax- 
ation eliminated the boundaries, reconstructing a single- 
crystal. This procedure indicated that the 40-atom cell 
is too small for the 0-terminated interfaces. The LDFT 
and shell-model relaxations for the 60-atom supercell pro- 
ceeded leaving the innermost atoms in a bulk-like en- 
vironment. Again, the adequacy of the axial supercell 
size was checked by relaxing an 80-atom supercell and 
by noticing the small diff'erence of 2 % between the inter- 
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face energies of the 60- and 80-atom supercells. The final 
interfacial configuration obtained with the two supercells 
showed no significant difference (smaller than 0.04 A in 
all the atomic positions within the supercell) 

B. Energetics and local atomic structure 

The supercells for the metastable grain-boundary 
structures identified in the previous section were then 
fully relaxed with respect to the inter-granular separa- 
tion (excess volume) and to all atomic positions. The 
interface energy was calculated as the difference between 
the total energy of the supercell, and the total energy 
of an equal number of AI2O3 formula units in the bulk 
phase, divided by the total grain-boundary area in the 
supercell. 

1. G(Al) interface 

First-principles LDFT calculations for the 60-atom su- 
percell predict a small axial excess GB expansion (2% 
of flrho, where arho=5.086 A cf. Ref. |l^) of the Al- 
terminated boundary, which lowers the interface energy 
from 1.92 J/m^ (corresponding to 0% expansion) to 1.88 
J/m^ (for 2% expansion). This interface energy is a rel- 
atively high value, compared to the ones of the rhombo- 
hedral and basal twins, namely 0.63 J/m^ (Rcf. |l^ and 
0.73 J/m2 (Ref. |ll|). 

The local atomic environment of this interface, which 
we anticipate is the most ordered one among the three 
metastable structures considered in this study, is consid- 
erably more disordered than the ones at the rhombohe- 
dral or basal twins. For example the atomic arrange- 
ments on opposite sides of the GB are not equivalent 
(Figure ^). In terms of Al-0 coordination, the cations in 
Al(l) are 5-fold coordinated while the ones in Al(l)* are 
6-fold coordinated. The bond lengths are summarized in 
Table ^ and compared to the ones for bulk alumina. 

In bulk alumina, each Al atom is coordinated by 6 
anions at two characteristic distances of 1.84 and 1.96 A. 
A similar splitting in the interfacial Al-0 distances was 
observed ip_.thc most stable of the rhombohedral twin 
structures,!!^ but not at the present S13 twin interface, 
where the bond lengths of the shell of first neighbours 
are considerably more spread. 

Similarly, there is a difference in coordination num- 
bers between the oxygen atoms in 0(2) and 0(2)*, which 
are 4- and 3-fold coordinated respectively. All the other 
oxygen atoms are 4-fold coordinated. In particular, the 
atomic environment around the innermost O atom is well 
bulk- like with two bond lengths around 1.82 A and two 
around 1.95 A (cf. bulk values of 1.84 and 1.96 A). 

The effect of atomic relaxations at the interface, com- 
pared to the bulk volume, is described by the differences 
in interlayer spacings shown in Table O. The relaxation 



involves mainly the first three atomic layers on both sides 
of the interface. As in the case of the bond lengths, there 
is a small difference between the atomic environment on 
opposite sides of the GB, which is also reflected in the 
change in interlayer distance (Table Note again that 
the central atoms are practically unaffected by the relax- 
ations induced by the boundary. These results suggest 
that the boundary effects involve mainly a slab of w 5 A 
about the interface plane. 

The free surface corresponding to this interface, de- 
noted as FS(Al), was obtained by artificially forming two 
grains with Al-terminating surfaces. These can be de- 
scribed as obtained by cutting the bicrystal structure of 
Figure || between layers Al(l) and Al(l),* and by sep- 
arating the grains until they do not interact any more. 
This is equivalent to consider only one of the two grains 
in the supercell and repeat it periodically along with 
a large lattice constant. The structure was then relaxed. 
The change in interlayer spacing due to relaxation is de- 
scribed in Table ||. The outermost Al layer contracts 
inward, while the outermost O layer relaxes outward, so 
that the final surface configuration is 0-terminated. The 
corresponding calculated surface energy is 2.77 J/m^, a 
relatively high value compared to the ones of the basal 
surface (1.94 J/m^ obtained by us with the present ap- 
proach, which was also used for the basal twins in Ref. 
|ll|; 1.95 J/m2 Ref. ||; 2.13 J/m^ Ref. ||) and the rhom- 
bohedral surface (1.98 J/m^ calculated with the present 
approach in Ref. pQ) . 

2. G(0) interface 

The axial relaxation of the G(0) interface, calculated 
by LDFT, yields that the minimum energy of 2.44 J/m^ 
corresponds to an axial excess GB expansion of 6% of 
flrho- The change in interlayer spacing due to atomic re- 
laxation for the 60-atom cell is shown in Table O. It sug- 
gests that the relaxation involves a slab of « 7 Acentered 
at the interface plane. In this case, the atomic environ- 
ments on opposite sides of the boundary are equivalent. 
The boundary atoms 0(0) are the most affected by the 
relaxation procedure, but also the inner 0(2) and 0(3) 
relax quite substantially (displacements of 4.6 % of Oi-ho). 

The local atomic environment is summarized in Ta- 
ble |. Despite the large atomic relaxations, the inter- 
facial oxygen 0(0) are in a bulk- like environment, be- 
ing 4-coordinated with two bond- lengths at 1.84 A and 
two at 1.99 A (cf. 1.84 A and 1.96 A in bulk). On 
the contrary, the next layers 0(1) and 0(2), and the in- 
ner layer 0(4) are 3-fold coordinated. (The fourth bond 
from 0(4), which points towards the interface, is elon- 
gated to 2.11 A. This is the case also for the 80-atom 
cell. Hence, the 60-atom cell is appropriate.) Similarly, 
the first three cations layers Al(l), Al(2), and Al(3) are 
under-coordinated with only 5 oxygen neighbors. In the 
80-atom supercell, the first atoms in bulk-like environ- 
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ment are the ones in Al(4), which are 6- fold coordinated, 
with three bond- lengths centered at 1.85 A and three at 
1.94 A. 

There are no obvious ways of defining a free surface 
from the interface configuration shown in Figure |. One 
possible variant, that has the advantage of creating stoi- 
chiometric and equivalent surfaces, is to assign half of the 
interfacial atoms 0(0) to one grain and half to the others. 
This creates a starting anion terminated configuration 
with a jagged structure. The structural minimization of 
this free surface leads to an oxygen terminated structure, 
indicated as FS(0), with a surface energy of 2.50 J/m^, 
hence being more stable than FS(Al). The changes in 
interlayer spacing due to relaxation are included in Ta- 
ble |l[ 

3. S(0) interface 

We already discussed that the initially imposed trans- 
lation state (0, 0.5) for this interface is not stable. Let 
us now specify this statement more precisely: the start- 
ing configuration (0, 0.5) identified with the shell-model 
calculations brings the oxygen atoms 0(0) and 0(1) very 
close together. We therefore started the LDFT structural 
optimization with a larger inter-granular distance T3 (6 
% of Oi-ho)- The atomic relaxation increased the distance 
between the "close" 0(0) and 0(1) layers, and this in- 
duced also a change in the lateral translation state. The 
minimum energy of 2.71 i/ir? is obtained at a rather 
large inter-granular separation, 12 % of aiho, which pro- 
duces a characteristic open boundary structure shown in 
Figure ^. In particular, there are large open channels 
along the boundary in the [2021] direction, with a typ- 
ical diameter of 6.7 A. Here impurity or dopant atoms 
of different species could in principle be accommodated 
interstitially without involving severe elastic strains. 

Atomic relaxations involve mainly the interfacial oxy- 
gen layer 0(0) and the first three atomic layers on both 
sides. The changes in interlayer spacings are summarized 
in Table |l[ These results suggest that the boundary af- 
fects mainly a slab of « 6 A centered at the interface. 

The innermost atoms Al(3) are practically unaffected 
by the interfacial relaxations, and thus they are in a bulk- 
like environment. The atomic distances of the first neigh- 
bouring O shell have the characteristic bulk-like splitting, 
three bond lengths being around 1.85 A and three around 
1.94 A. On the interface, cations and anions are respec- 
tively 5-fold and 3-fold coordinated. However, the anions 
of the second layer 0(2) have already the bulk coordina- 
tion number of 4. The Al-0 bond lengths are summarized 
in Table | 

The difference in interlayer separation between 0(0)- 
0(1) and O(0)-O(l)*, respectively 1.6 and 0.78 A(-t-8.5 % 
and -9.8 % of aiho, respectively, cf. Table |), suggests the 
O(0)-O(l) plane as an alternative cleavage plane to the 
ones already mentioned. The two free surfaces obtained 



by cleaving the bicrystal of Figure at these planes are 
not stoichiometric and geometrically not equivalent, but 
the supercell is overall stoichiometric. Both these FS con- 
tribute to an average surface energy of 2.76 J/m^. This 
average energy value provides some information about 
the possible cleavage mechanisms of alumina bicrystals, 
which we will discuss in the next section. Incidentally its 
value is close to the one of the surface energy of the sto- 
ichiometric Al-terminated surface (see section IV.B.l). 

C. Large inter-granular separations: bonding and 
cleavage 

Two sets of LDFT calculations monitored the interface 
energy as a function of large inter-granular separation in 
two different regimes. The starting points of the first 
one were the isolated grains with free surfaces, which 
were then brought closer, down to the formation of the 
three metastable interfaces previously described. These 
calculations are conceptually similar to the experimental 
procedure of diffusion bonding of bicrystals in ultrahigh 
vacuum, and should shed light on the boundary forma- 
tion mechanisms. The starting points of the second set 
of calculations were the relaxed metastable translation 
states of the bicrystals, described in the previous sec- 
tion, which were then gradually pulled apart, up to the 
formation of the free surfaces. These calculations reveal 
insights into the cleavage mechanisms. 

The fundamental difference between the two methods 
is that in the first case the cleavage plane is defined by an 
a priori choice of the free-surface terminating layer, while 
in the second case the cleavage plane is determined by the 
mechanical properties of the bicrystal, and therefore by 
the weakest bonds. 

The LDFT results for the interface energies versus the 
axial component T3 are shown in Figure |. The energies 
corresponding to the first set of calculations (from free 
surfaces to internal boundaries) are connected with solid 
lines, while the values for the second set (from internal 
boundaries to free surfaces) are connected with dashed 
lines. Solid horizontal lines represent the limit of free 
surfaces, being twice the energy of the free stoichiometric 
surfaces FS(Al) and FS(0), 5.53 and 4.94 J/m^ respec- 
tively. Each energy value corresponds to a fully relaxed 
configuration of atoms in the supercell. 

Let us first concentrate on intergranular bonding, i.e., 
the formation of the two most stable structures, G(A1) 
and G(0), from the corresponding free stoichiometric 
surfaces FS(Al) and FS(0), by reducing the intergranular 
distance. Note that the relative energetic stability of the 
free surfaces is reversed with respect to that one of the 
bicrystals: FS(0) is more stable than FS(Al), and G(A1) 
is more stable than G(0). With decreasing distance, the 
interface energies remain close to the values of the corre- 
sponding free surfaces down to an intergranular distance 
of 2.54 A (0.5 % of arho)- There a significant difference 
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of more than 5 % appears due to the onset of interaction 
between the two slabs, and from there the values mono- 
tonically decrease down to the equilibrium metastable 
configurations G(A1) and G(0). The crossover point in 
the relative energetic stability, when the Al termination 
becomes more stable than the O termination, is around 
1.8 A (0.35 % of 

For intergranular cleavage, starting from the equi- 
librium interfaces G(A1) and G(0), and stretching the 
bicrystals along T3, the reverse energy paths as for the 
bonding procedure are followed up to displacements of 
0.81 A (0.16 % of arho)- Subsequently, both the bicrystals 
react to the increase in intergranular distance by stretch- 
ing the weakest bonds: the cleavage plane passes through 
0(l)*-0(2)* in G(A1) and through O(0)-O(l) in G(0). 
Cleavage on these planes create non-stoichiometric free 
surfaces, leaving two terminating oxygen atoms on one 
side and one on the other side. 

At larger intergranular distances, the energies of the 
three metastable structures increase at the same rate, ap- 
parently to the same surface-energy limit. This may be 
understood by noticing that, in the limit of large sepa- 
ration, the two grains resulting from this 'asymmetric' 
cleavage in the G(A1) and G(0) supercells are struc- 
turally equal to the grains resulting from the cleavage 
through O(0)-O(l) in the S(0) supercell. Therefore, 
the corresponding three cleaved free surfaces become the 
same for large separation in the supercells. Their surface 
energies become equal to 2.26 J/rn? , which is the average 
surface energy of the two non-stoichiometric terminating 
surfaces identified for S(0) (see section IV. B. 3). 

D. Electronic structure 

The self-consistent LDFT results obtained for the three 
equilibrium metastable interface configurations were used 
to investigate the effect of the interfacial local atomic 
arrangements on the electronic structures of alumina 
bicrystals. In particular, we focus on the site-projected 
density of states (PDOS). We define the PDOS as the 
projection of the total density of states (DOS) on partial 
waves inside spheres centered on atomic sites. The radius 
of the spheres were 1.27 A and 0.58 A for oxygen and alu- 
minium sites, respectively. The crystal Bloch eigenstates 
were calculated on a 6 x 6 x 2 k-point mesh and broad- 
ened with Gaussians of 0.8 eV width. The same setup 
was used in the analysis . of . t |he rhombohedral E7 and 
basal E3 twin boundaries.t2Hlil 

The results are shown in Figure ^ where the oxygen- 
PDOS (0-PDOS) of the three interface states are com- 
pared with the 0-PDOS of bulk alumina. The latter is 
indicated by a thick solid line, while a dashed line and a 
thin solid line indicate the PDOS of oxygen atoms, close 
and far from the interface, respectively. To facilitate the 
discussion of the results, the interfaces PDOS have been 
rigidly shifted so that the bulk-like 0-PDOS (thin solid 



lines) is on top of the 0-PDOS for bulk alumina (thick 
solid line). The zero of energy in the plot is the Fermi 
level of the bulk supercell (the energy eigenvalue of the 
highest occupied valence state), and the Fermi levels of 
the supercells containing the interfaces are marked with 
vertical solid lines. 

A discussion of the Al-PDOS is not given here because, 
as already considered in Ref. it is less illustrative than 
the 0-PDOS in a-Al203: Due to the high ionic character 
of the Al-0 bond in alumina, all the features of the DOS 
(up to « 20 eV) are dominated by the 0-PDOS, the 
contribution of the Al-PDOS being 5—10 times smaller, 
and similar in the spectral shape of peaks and shoulders, 
even in the lower conduction band region. 

Let us first emphasize the similarity between the O- 
PDOS in bulk alumina (thick solid lines) with the ones in 
bulk-like environment far from the interfaces (thin solid 
lines). There are hardly any differences in the lower va- 
lence band (with energies centered around -17 eV), and 
minor differences in the higher valence band (with ener- 
gies centered around -5 eV) and in the conduction band. 
This shows that the effect of the local environment on the 
electronic structure is localized within the first 2-3 inter- 
facial oxygen layers, the most affected being always the 
interfacial ions (dashed lines) , as intuitively expected and 
already observed at the rhombohedral E7 twin boundary 
in Ref. |l^. This further confirms that the supercell di- 
mension along T3 is sufficiently large to leave the inner 
atoms in a bulk-like environment. 

We now discuss the differences between the interfa- 
cial and bulk-like valence 0-PDOS. Let us first notice 
that the E13 GB affect both the lower and higher va- 
lence bands of the interfacial O atoms. The differences 
are moderate in G(A1), where the interfacial 0-PDOS 
corresponding to 0(1) and 0(1)* (dashed lines in Fig- 
ure ^), have still the characteristic bulk-like structure 
(thin solid line). The biggest deviation appears for 
0(1),* which is 5-fold coordinated. More pronounced 
differences are present in G(0), where the interfacial O- 
PDOS (0(0) and 0(1)^0(1)*, dashed lines in Figure |b) 
are qualitatively different from the bulk-like structure 
(thin solid fine): the states corresponding to 0(1) (3- 
fold coordinated) are redistributed towards higher ener- 
gies, while the states corresponding to the 0(0) (4-fold 
coordinated) higher valence bands are redistributed to- 
wards lower energies, forming a two-peak structure. The 
high-energy S(0) boundary induces considerably more 
drastic changes in the electronic structure. The lower va- 
lence band is rigidly shifted by 1.8 eV to higher energies, 
and a three-peak spectral shape appears in the higher 
band. Let us also note that the valence-band edges for 
the G(A1), G(0), and S(0) supercells, respectively, are 
0.2 eV, 0.9 eV, and 1.3 eV higher than for bulk alumina. 

The changes in the unoccupied region of the 0-PDOS 
are now addressed. The main feature, common to the 
three interfaces, is the presence of a pre-peak in the en- 
ergy gap, at 6-7 eV, just below the conduction-band edge. 
A similar gap state was observed in the highest energy 
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structure of the rhombohedral twin. The fundamental 
band gap of bulk alumina (6.6 eV with the LDFT method 
used in the present calculations) is therefore reduced by 
the presence of this pre-peak to 5.41 eV in G(A1), 4.47 
eV in G(0) and 3.72 eV in S(0). 

The real-space projection along [1210] of the density 
of unoccupied electron states in the energy range be- 
tween 5 and 7 eV is shown in Figure ^ for the G(A1) 
interface, taken as the representative case. The density 
is proportional to the grey-scale, ranging from white for 
zero density, to black for the highest density. The pro- 
jected atomic positions in the supercell are marked by the 
labels Al and O. Note that most density of the electron 
states corresponding to the pre-peak is spatially confined 
at the interface, and that a considerable density fraction 
is localized on the first three interfacial O atoms. 

At higher energies, between 10 and 20 eV, the three 
interfaces modify the interfacial 0-PDOS substantially 
with respect to the bulk three-peak signal. Moreover, 
within the same boundary structure, there are signifi- 
cant differences between the 0-PDOS corresponding to 
different interfacial sites, like 0(1) and 0(1)* in G(A1), 
or 0(0) and 0(1) in G(0). This is likely to play an im- 
portant role in the interpretation of experimental ELNES 
spectra (cf, e.g. Rcf p^. 

As demonstrated in Rcf. |l3| for the rhombohedral S7 
twin boundary, the ELNES of the oxygen core-level ion- 
isation K edge, which is approximately proportional 
to the conduction band PDOS of angular-momentum 
p character at oxygen sites, can provide experimental 
informations to distinguish the local electronic states 
and atomic coordination at oxygen sites located in the 
bulk volume or at interfaces. The calculated 0-PDOS 
for the unoccupied conduction-band states of the E13 
twin boundary indicate characteristic differences between 
the three metastable structural models obtained by the 
LDFT structure optimisation. Hence, an ELNES exper- 
iment with a El 3 twin bicrystal in a scanning transmis- 
sion electron microscope with sufficiently high spatial and 
energetic resolutions, comparably to the ELNES experi- 
ment of Ref. |l^, appears very promising to augment our 
theoretical characterisation of this interface. 



V. CONCLUSIONS 

The energetic, structural, and electronic properties 
of the S13 (1014) mirror twin boundary in a-Al203 
were studied by atomistic modelling with both empirical 
shell-model potentials and first-principles LDFT calcu- 
lations. Three metastable variants of the interface were 
identified, one Al-terminated boundary, G(A1), and two 
0-terminated boundaries, G(0) and S(0). The lowest- 
energy interface is Al-terminated, as prpjiously predicted 
by empirical shell-model potentials .tj'E3 The interfacial 
energies of the fully relaxed interfaces, as calculated with 
LDFT, were discussed. 



Despite the thermodynamic stability of the G(A1) in- 
terface over the 0-terminating variants, a high-resolution 
transmission electron microscopy (HRTEM) experiment 
with high-purity diffusion-bonded bicrystals-detected the 
G(0) boundary, rather then the G(A1) oneJ23 In order to 
clarify these observations, the formation of the G(A1) and 
G(0) boundaries was studied by gradually reducing the 
intergranular distance, starting from the corresponding 
free surfaces FS(Al) and FS(0), to model theoretically 
the procedure of diffusive bonding in ultrahigh vacuum. 
This showed that the relative energetic stability between 
these two interfaces depends on the intergranular separa- 
tion, and is actually reversed for the free surfaces, FS(0) 
being more stable than FS(Al). From this a possible ex- 
planation for the HRTEM observations emerges: prior 
to bonding the free surfaces are 0-terminated. When 
the crystals are diffusion bonded, the equilibrium ther- 
modynamics would drive the atomic reconstruction of 
the 0-terminated interface into the Al-terminated one, 
G(A1). However, the energy barrier for such a recon- 
structive transformation is apparently too high, as long 
as no additional Al is provided in the bonding process, 
and the bicrystal is therefore locked in the metastable 
G(0) configuration. 

The analysis of the site-projected densities of states 
showed that the effect of the interfaces on the electronic 
structure is localized within 5-7 A around the boundary 
plane, and affects mostly the upper valence band and the 
conduction bands. In particular, the three variants of the 
S13 boundary narrow the band gap by defect states be- 
low the conduction-band edge, which are highly localized 
at the interface. 

An important final conclusion of this work is that, al- 
though the S13 twin boundary has a short lateral period- 
icity (which enables the thorough treatment by LDFT), 
the resulting structural complexities, as well as the inter- 
face energies of its metastable configurations, are signif- 
icantly higher than those of the previously investigated 
rhombohedral S7 and basal E3 twin boundaries, usu- 
ally called "special" grain boundaries. The E13 interface 
is therefore a promising model case for "general" grain 
boundaries in a-Al2 03. We anticipate that it will be 
useful to provide a theoretical understanding of the in- 
fluence of segregated cation impurities on the structure 
and energetics of alumina interfaces. 
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TABLE I. Equilibrium Al-0 distances in the interfacial 
and bulk-like regions of the G(A1), G(0), and S(0) super- 

ccUs. For comparison, there arc two Al-O distances in bulk 
a-Al203, with calculated values of 1.84 A and 1.96 A. 
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FIG. 1. Corundum structure of a-Al2 03 viewed from the 
[1210] direction. Oxygen and aluminium atoms are repre- 
sented by light-gray and black circles, respectively. Dashed 
lines mark the (1014) planes leading to stoichiometric S13 
interfaces with Al (B) and O (C) terminations, respectively. 




FIG. 2. Shell-model total energy (Ry/supercell) for the 
cation terminated interfaces as a function of the lateral trans- 
lation state Tiei -|-T2e2, with ei=[2021] and e2=[1210]. 
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FIG. 3. LDFT result for the relaxed atomic structure of 
the G(A1) interface viewed along the [1210] direction. The 
atoms are represented as in Figure |l|, and the labels on the 
right-hand side identify the atomic layers. The dashed box 
indicates the axial size of the 60-atom supercell in the [5054] 
direction, and both the size and lateral translation-state in 
the [2021] direction. The supercell contains two equivalent 
horizontal interfaces, one in the middle (marked by the dotted 
line), and one both at the top and bottom ends (marked by 
the dashed line segments) arising from the periodic boundary 
conditions in the [5054] direction. 
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FIG. 4. Shell-model total energy for the anion termi- 
nated interfaces as a function of the lateral translation state 
Tiei +T2e2, where ei = [202i] and e2 = [1210]. 
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FIG. 5. LDFT result for the relaxed atomic structure of 
the G(0) interface viewed along the [1210] direction. Atoms, 
labels, and symbols as in Figure ^. 
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FIG. 6. LDFT result for the relaxed atomic structure of 
the S(0) interface viewed along the [1210] direction. Atoms, 
labels, and symbols as in Figure ^. 
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FIG. 7. Interface energies (in J/m^) versus the axial inter- FIG. 8. a) 0-PDOS for the G(A1), b) for the G(0), and 

granular separation (in % of Orho) for the three metastable c) for the S(0) interfaces. The zero of energy is the Fermi 

interfaces G(A1), G(0), and S(0). The solid (dashed) lines level (energy eigenvalue of the highest occupied valence-band 

connect the results obtained by reducing (increasing) the in- state) of the bulk superceU, the interfaces PDOS have been 

tergranular separation from free surfaces (internal interfaces) rigidly shifted so that the bulk-like O-PDOS (thin sohd hnes) 

down to the the formation of the internal boundaries (free sur- is on top of the O-PDOS for bulk alumina (thick solid line), 

faces). Solid horizontal lines represent twice the energies of The Fermi levels of the supercells containing the interfaces 

the free surfaces FS(Al) and FS(0) (cf. Table 0; NB: upon are marked by solid vertical lines, 
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surface energy because the resulting surface area is twice the 
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FIG. 9. Real-space projection along the [1210] direction of 
the density distribution of unoccupied electron states in the 

energy range between 5 and 7 eV for the G(A1) interface. 
Al and O labels mark the projected atomic positions in the 
supercell. 
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